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I. INTRODUCTION 

The recent discovery of the Higgs boson at the LHC 
has completed the list of the known particles, providing 
the last missing element necessary for the standard model 
(SM) to be the framework for particle physics. However, 
it is well known that the SM suffers from various inter¬ 
nal problems indicating that this is not a fundamental 
theory, and in fact it should be considered just as an ef¬ 
fective low-energy theory. The strong CP problem is one 
of these problems. It emerges from adding to the QCD 
Lagrangian the so called 0 term 

£D0gTr(GAG), (1) 

written in terms of the QCD gluon field strength 2-form 
G. This is a renormalizable and gauge invariant term, 
which violates CP and is allowed in any generic gauge 
theory in four dimensions. In the SM it contributes to 
GP-odd observables such as the neutron electric dipole 
moment, which is stringently constrained by experiment, 
pushing the 9 parameter down to 10“^°. Since the nat¬ 
ural value of this parameter should be of order one, this 
becomes a fine-tuning problem. The question of why it 
turns out to be so small is the strong CP-problem. 

A solution of the strong CP problem has been found 
by Peccei and Quinn in the periodicity of the nonpertur- 
bative QCD 6 vacuum [1] by promoting the 6 parame¬ 
ter in Eq. (1) to be a held 9{x). Then the interaction 
9{x) Tr (GaG) generates in the 9 vacuum a nontrivial 
potential for 9{x), selecting a zero vacuum expectation 
value {9) = 0. The huctuations around this vacuum rep¬ 
resent a pseudoscalar held a(x), dubbed the axion. Then 
dynamically the GF-violating term (1) is replaced by the 
GP-conserving interaction a(x) Tr (GaG). 


The 9 parameter can be promoted to be a held, by 
means of a pseudoscalar held, (/>(x), of any origin, cou¬ 
pled to the Pontryagin density Tr (GaG) of the gluon 
held. This could be a Goldstone boson of a C(1)a sym¬ 
metry, spontaneously broken at some scale much larger 
than the electroweak scale of 250 GeV, to be compatible 
with the experimental data as well as with astrophysics 
and cosmology. There are many symmetry based propos¬ 
als of this kind in the literature, as possible solutions of 
the strong CP problem (for a recent review see Ref. [2]). 
A characteristic feature of this approach is that all the 
couplings of the axion are determined by the scale of 
symmetry breaking, which is a free parameter. 

On the other hand it is well known that various scenar¬ 
ios for the Planckian physics involve axionlike fields [3-6] . 
Those fields can play the same role as the conventional 
Goldstone-type axions in the solution of the strong CP 
problem, but with all their couplings completely deter¬ 
mined by the Planck scale. 

In particular the axionlike fields may appear rather 
naturally in a field theory on the torsionful manifolds 
with its metric sector treated as a “rigid” background. 
The first scenario of this kind was proposed in Ref. [7], 
where an axionlike field appears as a consequence of the 
constraint imposed on the quantum theory requiring the 
conservation of the torsion charge, as suggested by the 
classical theory. 

Recently, in Ref. [8], the axion has been introduced as 
a pseudoscalar field, the so-called Barbero-Immirzi (HI) 
axion, interacting with gravity via the Nieh-Yan (NY) 
density [9, 10]. One of the motivations for the intro¬ 
duction of this field was the possibility of eliminating 
the confusing divergence present in the C(1)a rotated 
fermion measure of the Euclidean path integral on the 
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manifolds with torsion. In addition to the usual Pon- 
tryagin density, in this case there appears a Nieh-Yan 
term, which becomes divergent when the regularization 
is removed [11]. The significance of this divergence was 
debated in the literature [12, 13] and a consensus on its 
status has not yet been reached. In the model of Ref. [ 8 ] 
this divergence can be absorbed by a redefinition of the 
Barbero-Immirzi field. This axionlike field was also pro¬ 
posed in Ref. [14] , in order to solve the strong CP problem 
in the Peccei-Quinn spirit. 

In the present paper, we show that the conservation of 
the torsion charge, within the framework of Ref. [7], is 
equivalent to demanding a vanishing Nieh-Yan density. 
This constraint can be implemented into the quantum 
theory by means of a Lagrange multiplier, identified with 
the so-called the Kalb-Ramond (KR) axion [15], due to 
its similarity with the axionlike held coming from string 
theory. 

Despite the starting points of Refs. [7]-[ 8 ] seeming to 
be different, they have the same physical properties when 
the torsion is integrated out. Therefore, within the ef¬ 
fective theory, the Kalb-Ramond [15] and the Barbero- 
Immirzi [ 8 ] axions are equivalent. We rigorously demon¬ 
strate this equivalence and study the solution of the 
strong CP problem based on these torsion-descended 
(TD) axions. Then we examine their possible cosmo¬ 
logical and astrophysical implications. 

In the present manuscript we concentrate on the dis¬ 
cussion of axions in Einstein-Cartan theory of gravity 
with torsion. For discussions on the role of axions moti¬ 
vated by Chern-Simons-type terms, see Ref. [16], where 
as a cosmological application, the accelerated expansion 
of the Universe has been considered [16, 17]. 

The article is organized as follows. In Sec. II we sketch 
the classical Einstein-Cartan gravity with torsion. In 
Sec. Ill we discuss the two procedures for quantizing the 
model, and show that despite their different origin, they 
are equivalent as effective theories. In Sec. IV we con¬ 
sider some cosmological and astrophysical implications 
of these torsion-descended axions. Finally, we summa¬ 
rize our present study in Sec. V. 

II. CLASSICAL GRAVITY SETUP 

We consider the Einstein-Cartan theory of gravity - 
a minimal construction of gravitational theory allowing 
the connection to possess a nonvanishing torsion [18, 19]. 
The Cartan’s structure equations, relating curvature and 
torsion with the vierbein e“ and the spin connection 
read 

de“+ u;“ 6 Ae''= r“, ( 2 ) 

da;“c + A (3) 

where e“ = e^dx^ and are the vier¬ 

bein and spin connection 1-form respectively. Hereon, 
bold symbols will denote differential forms, while greek 


and latin indices stand for spacetime and Lorentz indices 
respectively. Additionally, the vierbein is related with 
the curved metric g^u through 

9p,u = gab e“ el , (4) 

where gab = diag (—1,4-1,-1-1,-1-1) is the Minkowski met¬ 
ric in four dimensions. The spin connection can be 
split into -|- K.°‘^, where the tilde indicates 

a torsion-free quantity and the contorsion tensor 
encodes the information about the torsion through the 
relation 7”“ = /C“{,Ae*’. Similarly, Eq. (3) together with 
the decomposition of the spin connection yields to 

^ac = ^ + /C“6a/C^" . (5) 

The Einstein-Cartan action can be written as 

5gr = ^ / Cabcd , (6) 

where = 8ttGn = SnMpi with Gn and Mpi being the 
Newton’s constant and the Planck scale, respectively. 

In the following, the SM fields are assumed to live in a 
curved torsionful spacetime. The nontrivial coupling of 
matter with torsion enters in the fermionic sector through 
the covariant derivative, 

T>iIj= dil}+^u:°-^^abtlJ + ^eA'^p + igBtjj, (7) 

where A and B denote the U (l)em and SU{3)c gauge bo¬ 
son 1-forms respectively and jab = 5 [ 7 a, 7 b]- Therefore, 
the complete model is described by the action 

>5 = 5gr - i ^ y {tpf 7A*X>-0/ - T>tpfA*-yiljf) 

~\j j Tr[GA*G] , ( 8 ) 

where '0 = —is the usual Dirac adjoint, 7 = 7 a e“ 
and the subscript / indicates the SM fermionic flavors. 
The symbol * denotes Hodge duality, while G and 
are the SU{3)c and U{l)em gauge field strength 2 -form 
respectively. 

When the action in Eq. ( 8 ) is varied with respect to the 
vierbein field, one obtains the corresponding Einstein- 
Cartan equation of motion, 

P'ab '^gabP- — ^ "^ab; (9) 

where Tab is the energy-momentum tensor of the system. 
In its form Eq. (9) looks similar to the Einstein equations 
derived within general relativity (GR). However, the two 
equations are different because the presence of torsion 
gives rise to an antisymmetric part in both sides of the 
equation, absent in GR. 

On the other hand, the variation of the action with re¬ 
spect to the spin connection yields to an algebraic equa¬ 
tion of motion for the torsion, 

Tabc = —, ( 10 ) 
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where = i'i/j/ 7 ‘^ 7 ^' 0 / is the fermionic axial current. 
Notice that only the completely antisymmetric part of 
the torsion couples to the fermionic fields. This part 
corresponds to its axial irreducible component, 

S = j,eabcdr^^^e^. ( 11 ) 


The quantum generating functional, with this condi¬ 
tion incorporated, takes the form 

Z = 5 (d*S) , (13) 

<p 


Now we rewrite Eq. ( 8 ), showing explicitly only those 
terms which depend on the torsion axial component 

S = So + ^ J (12) 

where Sq = Sgr +S^ + 5gk and 5gk represents the gauge 
field kinetic terms [the last two terms in Eq. ( 8 )]. The 
entity = X)/('^/)o denotes the axial current 1 -form, 
where the sum runs for all the fermionic flavors /. 


III. QUANTUM THEORY AND TD AXIONS 


where 5 [(/ 5 ] is the action given in Eq. (12) and (p denotes 
all the fields except for e“, treated as a rigid background. 
The argument of the delta in Eq. (13) can be passed to 
the effective action using the integral representation 


6 (d*S) 


^cj) 


^ (14) 


This allows us to write 


Z = y J SiSSicj) exp I* 




Sa*S 


- ^J5a*S - d(j)AirS 


(15) 


The quantization of the model with the action given in 
Eq. (12) can be carried out on the basis of the path inte¬ 
gral representation for the generating functional. How¬ 
ever, at present it is unknown if this procedure is ap¬ 
plicable to the quantization of the whole gravity sector 
(for the current status of this problem see, for instance. 
Refs. [20]). 

In the scenario studied here, the SM fields lie on a 
torsionful manifold, whose only quantum gravity effects 
enter through the torsion, while the metric or Riemma- 
nian curvature remains as classical variable. Quantum 
torsion seems to be easily treatable because the equa¬ 
tion of motion for the torsion ( 10 ) is algebraic, showing 
that the torsion is a nonpropagating field, which can be 
exactly integrated out from the theory. 

However, as it was observed in Ref. [7], this treatment 
of the torsion should be done with caution. It follows 
from Eqs.(I0)-(II), that S oc J^. Since the action (12) is 
U{1)a symmetric, the Nother current Jj is conserved at 
classical level, leading, as follows from the above relation, 
to the conservation of the torsion charge Qs = / *S. 

On the other hand, we know that the fermionic mea¬ 
sure of the path integral is not U (1)^ invariant. This fact 
manifests as the anomalous nonconservation of J® at the 
quantum level. As pointed out in Ref. [7], this must be 
taken into account before integrating out the torsion, in 
order to maintain the self consistency of the constructed 
effective theory. Following Ref. [7], an effective quantum 
theory can be constructed through a constraint requiring 
the conservation of the torsion charge d*S = 0. Notice 
that this is a gauge invariant condition, which is impor¬ 
tant for the self consistency of the SM sector of the the¬ 
ory. Later we show that this condition eliminates the 
divergent part of the U{\)a anomaly mentioned in Sec. 
I and affecting the tractability of the quantum theory in 
the presence of the torsion. 


where J® = J^. 

Since S is a nonpropagating field, we can integrate it 
out in the standard way [7, 19] carrying out a variable 
transformation 

S' = S--«;2d((.-iK2j5^ (jg) 

3 2 

with the Jacobian equal to unity. This new variable ap¬ 
pears in the exponent in Eq. (15) only in the bilinear 
combination Sa *S and, therefore, can be exactly inte¬ 
grated out. As a result we get the effective action 

5eff = 5o- ^ J 3^a*3^- 

-\j dchA^dch-f J Ci>d*J^ (17) 

For convenience we have made a redefinition $ = 
(j). Notice that the integration out of the torsion 
makes $(x) a dynamical field with the canonical kinetic 
term. As follows from the last term in Eq. (17) this field 
is pseudoscalar. It is what we called in the introduction 
KR axion field. 

At the quantum level, the last term of Eq. (17) is noth¬ 
ing but the axial anomaly [21]. In the path integral lan¬ 
guage d*J® 0 is the manifestation of the U{1)a non¬ 
invariance of the fermionic measure [ 22 ], mentioned in 
the Introduction. 

On the Riemann-Cartan manifolds the axial anomaly 
was first studied in Refs. [23]. However, as shown in 
Ref. [11], the computation of such an anomaly gives rise, 
in general, to an additional previously missed term, called 
Nieh-Yan topological density [9] so that under a U{1)a 
rotation of the fermion fields the fermion measure ex- 
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periences a nontrivial variation [11, 23] 

- r\J^ 

TT 

+ ^Tr[GAG] + ^7^“^7^a^, (18) 

+ 2M2 (raAr“ - e^AOb 

Here aem and are the electromagnetic and QCD cou¬ 
plings, respectively, Nj is the total number of fermionic 
flavors, Nq is the number of quarks and 
where Qf is the charge of / fermionic flavor. The last 
term in Eq. (18) is NY topological density with the regu¬ 
lator multiplier being divergent when the regularization 
is removed, M —>■ oo. As mentioned in Sec. I the status 
of this divergence is still debated in the literature. 

However we find that in the approach of Ref. [7] it is 
irrelevant since the NY term AT vanishes identically due 
to the condition d*S = 0 imposed on quantum theory by 
insertion of the corresponding delta function in Eq. (13). 
In fact this follows from the identity derived in Ref. [9], 

at = r“ A Ta - Tlab A e“ A e*- = d(e“ A Ta ), (19) 

and the definition of the field S in Eq. (11) written in 
the form 



S)Ip Slip Ip Slip X exp 


ha I 


In the effective action (23) we added the QCD and the 
gravitational 6 and 0 terms. They are the gauge and 
gravitational Pontryagin densities allowed by the gauge 
symmetries of the theory. These terms are also needed 
for the model completeness, and play the role of countert¬ 
erms for the axial anomaly quantum corrections. They 
do not affect the previous derivation, since due to their 
topological nature they do not change the equations of 
motion. 

Recently in Ref. [8] there has been proposed an alter¬ 
native scenario in gravity with torsion also leading to an 
axionlike field. This scenario is inspired, in particular, by 
the Chern-Simons modified gravity motivated in its turn 
by string theory. The gravitational action according to 
Ref. [8] is modified at the classical level by the term 

Stot=S +J Pix)Ar, (25) 

where S is the action given in Eq. (8) or (12). This action, 
being used in the quantum generating functional, allows 
one to absorb the divergent NY part of the anomalous 
17(1)A variation of the fermion measure (18) by a renor¬ 
malization of the field pi{x) called in Ref. [8] the HI axion. 
The field pi{x) becomes a dynamical field with the canon¬ 
ical kinetic term after excluding a nondynamical torsion 
held using the classical equation of motion 


*Soce„Ar“. (20) 

Then from Eqs. (19)-(20) it follows that 

dxS = 0 A/" = 0. (21) 

Thus, neglecting the Nieh-Yan term in the axial anomaly, 
we can write for the axial current 


d * J5 = _ xr [GaG] 


27T 


Nf ~ab 


( 22 ) 


87r2 


TZ /\TZab- 


The right-hand side is written in terms of torsion-free 
quantities. This is attainable by the introduction of 
proper counterterms, as shown in Ref. [7]. 

Now we substitute the identity (22) into Eq. (17) and 
obtain the resulting effective action of the model. 


Here we introduced a parameter 

/<(, = K-i - 4 X 10^® GeV, (24) 


analogous to the decay constants of helds with derivative 
couplings, such as Goldstones of spontaneously broken 
symmetries. 


S=^«,2d^+1^2j5^ (26) 

derived from the action (25). This is equivalent to in¬ 
tegrating out the torsion held S in the generating func¬ 
tional carrying out the transformation (16). Note that 
the held /3(x) in the classical action (25) is nothing but a 
Lagrange multiplier setting the classical level constraint 
A7 = 0. Now one can immediately realize that in view 
of the identities (19))-(20) it is equivalent to the con¬ 
straint d*S = 0, which in the approach of Ref. [7] was 
set at quantum level as a constraint incorporated in the 
generating functional (13) or (15). On the other hand 
both approaches lead to the same effective quantum the¬ 
ory with the effective action given in Eq. (23) with the 
identihcation of the KR and HI axions <l>(a;) = /3(a:). As 
we have seen KR and HI axions originate from rather 
different treatments of quantum theory in the presence 
of the torsionful gravity. Nevertheless from the point of 
view of low-energy effective theory and the resulting phe¬ 
nomenology they both are equivalent particles, which we 
call from now on TD axions. 

Additionally, the TD axions may also appear in the 
context of the torsion-induced quintessential axions [16]. 
In this framework, the axial current is modihed by the 
addition of the Chern-Simons-type terms, in order to be 
conserved in the zero mass limit. The complete cancella¬ 
tion of the torsion sector in the anomaly can be addressed 
requiring the torsion to be an exterior derivative of a 
pseudoscalar held, identihed later with the axion [16]. 
Remarkably, this approach leads to the same effective 
theory as Refs. [7, 8] for a constant dilaton held. 
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IV. PHENOMENOLOGY AND COSMOLOGY 
WITH TD AXIONS 

In the effective action in Eq. (23) of the considered TD 
axions the last term is the most important in the context 
of the strong CP problem. The presence of the coupling 
of an axionlike field to the gluon field Pontryagin density 
is the necessary and sufficient condition for the solution of 
the strong CP problem via the Peccei-Quinn mechanism. 
The TD axion decay constant /$ introduced in Eq. (23) 
represents a typical energy scale of the model, related to 
the Planck scale. Thus, the TD axions $ emerge without 
any accompanying free parameter. This drastically dis¬ 
tinguishes them from the axions introduced as Goldstone 
fields of spontaneously broken symmetries, requiring at 
least one free model parameter, i.e., the scale of symme¬ 
try breaking. 

In principle in certain models both the TD and Gold- 
stone axions can coexist mixing with each other [14] . We 
do not consider this case here since in the presence of the 
TD axions, solving the strong CP problem without free 
parameters, introduction of other axions looks excessive. 

Focusing on the last term in Eq. (23), let us write down 
9{x) = 9 + $(x)iVq//$. The main point of the Peccei- 
Quinn mechanism is that the coupling ~ f 9{x) Tr [GaG] 
generates a nontrivial potential for the 9{x) field. The pe¬ 
riodicity of this potential in 9 [24] selects the unique non¬ 
trivial minimum {9{x)) = 0, corresponding to (<I>(x)) = 
—9f^/Nq. Perturbations around this vacuum generate 
the physical pseudoscalar axion field a{x) = $(x) — 
(<i)(x)). Thus the only surviving piece of the last term of 
Eq. (23) is the GP-conserving interaction a{x) Tr [GaG]. 
This solves the strong CP problem with the help of the 
TD axions. The mass of the TD axion can be calcu¬ 
lated in the usual way, as is done for any axion field (for a 
review cf. Ref. [25]), and depends only on the parameter 
f<s> defined in Eq. (24). The nontrivial mass is generated 
by instantons and for the value in Eq. (24) it turns out 
to be 


ma 


nriT, 


U ^muTUd 
U 'mu + nT-d 


10-12 eV, 


(27) 


where /tt = 93 MeV is the 7 r-meson decay constant and 
mT^,mu,d are the masses of 7r-meson and u,d quarks. 
Such an extremely light particle, having the inverse 
Planck mass suppressed interactions with gauge fields, 
is unobservable in laboratory experiments. 

Nonetheless, an axion with these properties may play 
a significant cosmological and astrophysical role, since 
it must satisfy the existing limits related to its origins. 
These aspects of the torsion-descended axions considered 
here, have been studied in Ref. [14]. It has been shown 
that such axions safely pass all the known astrophysical 
constraints, which originate from the energy loss of a 
stellar core in the form of axions and its impact on stellar 
evolution. 

An interesting cosmological prediction applied to the 
torsion-descended axions, also discussed in Ref. [14], 


is the production of axion isocurvature perturbations 
in which amplitude is constrained from the above by 
WMAP data [26]. Assuming that they are the domi¬ 
nant component of dark matter in the Universe it was 
found for the upper limit [14] 

i7/<10i°GeV, (28) 


for the Hubble expansion rate Hj during inflation. Now 
we can estimate the tensor-to-scalar ratio r = Pt/Ps 
using expressions for the power spectra of the scalar Ps 
and tensor Pt perturbations [27] 


Ps 


1 




(29) 


Here Mpi = I/^/SttGn ~ 2.44 x 10^® GeV is the reduced 
Planck scale and e is the standard slow-roll parameter. 
We use the Planck Collaboration result Ps ~ 2.19 x 10“® 
(see Refs. [28, 29]) and find from Eqs. (29) 

Hi ft! 2.5 X 10i'‘ GeVyf. (30) 


Using the limit (28) one finds the prediction for the cos¬ 
mology with torsion-descended axions 

r< 1.6x10-®. (31) 


This is in dramatic contradiction with the recently pub¬ 
lished result by BICEP2 [30], r = 0.2^q q 5 . Nevertheless, 
the situation has recently changed after the publication 
of the Planck Collaboration detailed analysis of the im¬ 
pact of the diffuse galactic dust polarized emission on 
the measurements of the polarization of the cosmic mi¬ 
crowave background (CMB) [31]. It has been shown that 
the BICEP2 result can be accounted for in the presence 
of this dust. Thus, the torsion-descended axions, con¬ 
sidered in the present paper, are currently not excluded 
by the cosmological data. Their cosmological test via 
the tensor-to-scalar ratio r is postponed for the future. 
Results of improved measurements of r, taking into ac¬ 
count the complications with the diffuse galactic dust, 
are expected to come in the near future from the Keck 
Array [32] and BICEPS [33] telescopes. The first results 
of the joint analysis of BICEP2/Keck Array and Planck 
data have been recently issued [34] showing only an up¬ 
per limit r < 0.12 at 95% C.L. New results from BICEPS, 
which will improve this limit, are expected during 2015 
and 2016 seasons [34]. 

The following final note might be in order. If the 
considered scenario is incorporated into the extradimen- 
sional setup [35] amended with the torsion [36] the fun¬ 
damental D-dimensional Planck scale M* could be re¬ 
duced down to the TeV values, dramatically changing 
the phenomenology and cosmology of the TD axions. 
In fact, making a rescaling Mpi —M* > 100 GeV in 
Eqs. (24) and (27), we find the values rua < 38keV and 
/$ > 100 GeV. Then for the rate of a —77 we get 
ra 77 < 10-^®eV. As to the cosmological aspects of the 
extradimensional TD axions, they require a dedicated 
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study. In particular, the value of the tensor-to-scalar ra¬ 
tio r cannot be obtained from (31) by the simple rescaling 
of Mpi. Let us recall that the bound in Eq. (31) was ob¬ 
tained using the limit (28) derived in Ref. [14] with the 
assumption that the axions be out of thermal equilibrium 
with photons during inflation. This condition can be vi¬ 
olated for the values of given above. The role of the 
extradimensional TD axions as a dark matter candidate 
should also be reconsidered. The corresponding study is 
in progress and its results will be published elsewhere. 

V. CONCLUSIONS 

We considered a solution of the strong CP problem 
via the Peccei-Quinn mechanism, implemented into the 
theory of gravity with torsion. We showed that the 
self-consistency condition of quantum theory d*S = 0 
proposed in Ref. [7] is equivalent to the requirement of 
vanishing Nieh-Yan topological density on the spacetime 
manifold. The Lagrange multiplier field, incorporating 
this constraint, leads to the torsion-descended axion cou¬ 
pled to the gluon Pontryagin density, Tr [GaG] , and 
therefore allows application of the Peccei-Quinn mech¬ 
anism for solving the strong CP problem. 

We considered the Kalb-Ramond and the Barbero- 
Immirzi axions proposed in the literature from quite 
different theoretical perspectives. We found that from 
the viewpoint of the effective theory these two torsion- 


descended axions are equivalent. 

An important property of the torsion-descended ax¬ 
ions is that their phenomenology has no free param¬ 
eters, rather they are completely determined by the 
Planck scale or, equivalently, by Newton’s gravity con¬ 
stant. The torsion-descended axion masses and their 
characteristic decay constants, are extremely small due 
to the Planck suppression, typical for this family of ax¬ 
ions rooted in gravity. We demonstrated the compati¬ 
bility of the torsion-descended axions with all the exist¬ 
ing cosmological and astrophysical limitations, as well as 
prospects for testing them in the near future measure¬ 
ments of the tensor-to-scalar ratio of the perturbation 
modes of the CMB. We also estimated the possible role 
of extra dimensions in phenomenology and cosmology of 
torsion-descended axions. 
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